A Model of Metallicity Evolution in the Early Universe 

G. J. Wasserburg^ and Y.-Z. Qian^ 

Received ; accepted 



^The Lunatic Asylum, Division of Geological and Planetary Sciences, California Institute 
of Technology, Pasadena, CA 91125. 

^School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455; 
qian@physics.umn.edu. 



- 2 - 



ABSTRACT 

We apply the phenomenological model used to explain the abundances of 
Fe and r-process elements in very metal-poor stars in the Galaxy to [Fe/H] 
of damped Lya systems. It is assumed that the first stars formed after the 
Big Bang were very massive and promptly enriched the interstellar medium 
to [Fe/H] ~ —3, at which metallicity formation of normal stars took over. 
Subsequent Fe enrichment was provided by Type II supernovae. The range of 
[Fe/H] at a given redshift z for damped Lya systems is explained by the time 
t* after the Big Bang at which normal star formation started in an individual 
protogalactic system. The average is ~ 80% the age of the universe for 
damped Lyo; systems at 2; ~ 1.5 to 4.5, indicating a long delay between the Big 
Bang and the turn-on of protogalaxies. It is inferred that a substantial fraction 
of the total baryonic matter may not have been aggregated into protogalaxies 
where normal star formation had occurred down to z ~ 1.5. The data near 
z — Ifl suggest that the rate of turn-on of protogalaxies was initially very low 
and slowly reached a maximum at ~ 3 Gyr after the Big Bang. This may be 
important in understanding the rate of formation of quasars. 

Subject headings: galaxies: abundances — galaxies: evolution — quasars: 
absorption lines 
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1. Introduction 

In a recent paper Prochaska & Wolfe (2000) reported new data on [Fe/H] for damped 
Lya systems at high redshifts. They also summarized available data with redshifts from 
^ fti 1.5 to 4.5 (e.g., Lu et al. 1996; Lu, Sargent, & Barlow 1997; Prochaska & Wolfe 1999). 
Their observations are in general accord with previous studies and expand the database at 
z > 3 substantially. Prochaska & Wolfe (2000) emphasized that there is a wide spread in 
[Fe/H] at a given z and no damped Lya system has [Fe/H] < —2.7 (see also Lu et al. 1996, 
1997). They also pointed out that there is only a minimal growth of [Fe/H] from z ~ 4.5 
to 1.5. In a theoretical study of Fe and r- process abundances in very metal-poor stars in 
the Galaxy by Wasserburg & Qian (2000, hereafter WQ), it was proposed that the first 
stars formed after the Big Bang were very massive 100 Mq) and promptly enriched the 
interstellar medium (ISM) to [Fe/H] ~ —3, at which metaUicity formation of normal stars 
(with masses ~ 1-60 Mq) took over. Subsequent Fe enrichment was provided by a subset of 
Type II supernovae. The interpretations of WQ were based on observations of metal-poor 
stars in the Galaxy by Gratton & Sneden (1994), McWiUiam et al. (1995), McWiUiam 
(1998), and Sneden et al. (1996, 1998). The apparent agreement between the lower bound 
on [Fe/H] of damped Lya systems and the critical metaUicity [Fe/H] ~ —3 deduced by WQ 
for transition from formation of very massive stars to normal stars suggests that this Fe 
enrichment model deserves further study. 

Here we present a phenomenological model for the range of [Fe/H] at a given z for 
damped Lya systems based on the model of WQ. It is assumed that formation of normal 
stars started in a damped Lya system (i.e., a protogalactic system turned on) at a time t* 
after the Big Bang. The assembly of [Fe/H] at a given z is interpreted as a sampling of 
t* ranging from to the age of the universe at z, t{z). The value ~ corresponds to 
the upper bound on [Fe/H]. Furthermore, the distribution of t* at a given z indicates the 
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rate of turn-on of protogalaxies prior to t{z). The data near z — 2.2 suggest that this rate 
was initially very low and slowly reached a maximum at ~ 3 Gyr after the Big Bang. We 
describe the Fe enrichment model in detail in §2 and apply it to explain the data on [Fe/H] 
of damped Lya systems in §3. Discussion and conclusions are given in §4. 

2. Fe Enrichment and Abundances in Metal- Poor Stcirs 

The Fe enrichment model of WQ was developed to explain the relation between 
abundances of Fe and r-process elements (r-elements) in metal-poor stars in the Galaxy. 
Meteoritic data on the inventory of radioactive ^^^I and ^^^Hf in the early solar system 
require at least two distinct Type II supernova sources for the r-process (Wasserburg, Busso, 
& Gallino 1996; Qian, Vogel, & Wasserburg 1998; Qian & Wasserburg 2000). These are 
the high-frequency H events responsible for heavy r-elements with mass numbers A > 130 
(e.g., Ba and Eu) and the low-frequency L events responsible for light r-elements with 
A < 130 (e.g., Ag). The recurrence timescales for the H and L events are Ah ~ 10^ yr 
and Al ~ 10^ yr, as required by replenishment of the appropriate radioactive nuclei in 
a standard mixing mass (~ the size of a molecular cloud) for a supernova. Additional 
evidence in support of different sources for the heavy and light r-elements has been found 
by Sneden et al. (2000). 

The observed wide dispersion in abundances of the heavy r-elements such as Ba and 

Eu over a narrow range of [Fe/H] 3 to -2.8 (McWiUiam et al. 1995; McWiUiam 1998; 

Sneden et al. 1996, 1998) led WQ to conclude that the H events cannot produce a significant 
amount of Fe. In contrast, there is a correlation between abundances of Fe and the heavy 
r-elements at [Fe/H] ;^ -2.5 (Gratton & Sneden 1994; see also McWiUiam et al. 1995). 
As Type la supernovae would occur only at metallicities much higher than [Fe/H] = —2.5, 
Fe enrichment of very metal-poor stars with [Fe/H] ^ —2.5 must be provided by the L 
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events. Consequently, the L events are responsible for the part of Fe contributed by Type 
II supernovae in general. There are ~ 100 L events during the time of ~ 10^° yr prior to 
solar system formation. To provide fti 1/3 of the solar Fe inventory by these events, each L 
event must enrich a standard mixing mass with [Fc/HJl ~ —2.5. 

The near absence of the heavy r-elements in stars with [Fe/H] ~ —4 to —3 (McWilliam 
et al. 1995; McWilliam 1998) and the sharp increase in the abundances of such elements at 
[Fe/H] ~ —3 to —2.8 led WQ to conclude that a source other than Type II supernovae must 
exist to produce Fe (along with other elements such as C, N, O, Mg, and Si) at [Fc/H] <^ —3. 
This source was attributed by WQ to very massive stars (with masses ^ 100 Mq) that first 
formed after the Big Bang. They further argued that formation of normal stars could not 
occur until [Fe/H] ~ —3 was reached. Presumably, this critical metallicity corresponds to 
conditions in the ISM that permit sufficient cooling and fragmentation to occur in collapsing 
gas clouds. A recent study by Bromm, Coppi, & Larson (1999) suggests that the very first 
stars were rather massive. The products of very massive stars formed from Big Bang debris 
have been discussed earlier by Ezcr & Cameron (1971). However, nucleosynthesis in such 
stars remains to be tested with adequate stellar models. 



3. [Fe/H] of Damped Lya Systems 

To discuss [Fe/H] of a damped Lya system at a given z, we assume the following 
history for its evolution: (1) at time ti after the Big Bang, matter consisting of Big Bang 
debris was isolated to form a system; (2) prompt enrichment by the first very massive stars 
ended at time t2, resulting in [Fe/H] ~ —3 in the ISM; and (3) formation of normal stars 
began at time t*, with the L events providing further Fe enrichment to an average ISM at 
regular intervals. The above times are related as ti < t2 < t*. As t* < t{z), it is necessary 
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to define the value of t{z) that is used. We take the redshift z to correspond to a time 

t{z)^lH^'n-'/\i + z)-y^ (1) 

after the Big Bang. We take the Hubble constant Hq = 65 km s~^ Mpc~^ and the matter 
contribution to the critical density = 0.3 [equation (|l|) gives essentially the exact result 
for t{z) at z > 1.5 for a flat universe with = 0.3 and a cosmological constant]. By our 
assumption, for t{z) > t2 the metallicity of a damped Lya system is at least [Fe/H] ~ —3. 
Prochaska & Wolfe (2000) reported that no system at 2 ~ 1.5 to 4.5 has [Fe/H] < —2.7. 
This means that the time required to provide prompt Fe enrichment must be less than 
t{z ^ 4.5) ^1.4 Gyr. The lower bound on [Fe/H] due to prompt enrichment is shown as a 
band between —3 and —2.8 in Figure 1 together with the data summarized in Prochaska & 
Wolfe (2000). 

In general, for t{z) > t* the metallicity of a specific damped Lya system at z is given 

by 

(Fe/H) = (Fe/H)p + {Fe/E)J-^^, (2) 

where the number ratio in round brackets (Fe/H) is related to the standard square bracket 
notation by [Fe/H] = log(Fe/H) — log(Fe/H)0, with (Fe/H)p corresponding to the prompt 
enrichment and (Fe/H)L to the contribution from a single L event to a standard reference 
mass of hydrogen. We take [Fe/H]p = —2.8, [Fe/H]L = —2.5, and Al = 10® yr, the same 
parameters used by WQ. As recognized by Lu et al. (1996), the scatter in [Fe/H] at a given 
z for damped Lya systems might result from their different formation histories. Equation 
(§) explicitly states that the range of [Fe/H] at z is caused by the different start times t* 
for normal star formation in individual systems (see Figure 2a). The values for [Fe/H] 
corresponding to t* = 1 and 1.5 Gyr are shown in Figure 1. Damped Ly« systems that just 
turned on at t* ^ t{z) would have [Fe/H] ^ [Fe/H]p. An upper bound on [Fe/H] exists as 
damped Lya systems that turned on at t* ~ would have the longest history of normal 
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star formation, and hence, the highest [Fe/H]. This bound is insensitive to the choice of 
[Fe/H]p. It can be seen from Figure 1 that almost all the data lie below the upper bound 
(usually well below this bound). Even the three exceptions are close to this bound (see §4). 

At a fixed z, equation (^ with equation (|I]) can be used to determine the start time t* 
for normal star formation in a damped Lya system from its [Fe/H]. In turn, a histogram 
of the number of systems at z within a given [Fe/H] interval determines the probability 
p (t*, t(2;)) dt* for normal star formation to start in the interval between and t* + dt* after 
the Big Bang. Knowing the probability distribution p{t*,t{z)) over < t* < t{z) at z, we 
expect that the probability distribution a,t z' > z can be obtained by discarding the part of 
p {t*, t{z)) for t* > t{z') and renormalizing the remaining part over < < t{z') (see Figure 
2b). For the case of a simple power-law distribution p (i:*, t(z)) = [{a + l)/t{z)][t* /t{z)]" . In 
this case, the average start time for normal star formation in damped Lya systems at z is 
{t*iz)) = [ia + l)/ia + 2Mz). 

The data in Figure 1 show a relatively high concentration in the interval 2.0 < z < 2.4. 
The average start time in this interval is {t*{z = 2.2)) ^ 2.5 Gyr. The frequency of 
occurrences of t* calculated from the data in this interval is shown as a histogram in Figure 
3. It can be seen that the frequency of occurrences is low for t* ~ (close to the Big Bang) 
and slowly increases to a maximum at t* ~ 3 Gyr. Assuming a power-law distribution for 
t*, we obtain a ^ 3, for which {t*{z)) ~ 0.8t{z) (see the corresponding curve for [Fe/H] 
shown in Figure 1). Values of t* for all the data are shown in Figure 4. The clustering of 
t* close to t{z) indicates that typically there is a long delay between the Big Bang and the 
start time t* for normal star formation in damped Lya systems. 



-8- 



4. Discussion and Conclusions 

We consider the dominant cause of dispersion of [Fe/H] at a given z for damped 
Lya systems to be the variation in the start time (t*) after the Big Bang for normal star 
formation in different protogalaxies. The bounds on [Fe/H] from our model appear to closely 
define the observed ranges. The average for damped Lya systems at z is ~ 0.8^(2;). 

The rate of turn-on of protogalaxies was initially very low and slowly increased to a 
maximum at ~ 3 Gyr after the Big Bang. It is not possible to identify a turnover in this 
rate without more data at 2; < 2. We suggest that the approach outlined here is a method 
for dating the start time for normal star formation in protogalaxies. As the formation of 
quasars is closely related to star formation, we consider that the histogram shown in Figure 
3 may offer some insights into the rate of formation of quasars. For example, if the rate of 
quasar formation it^g is proportional to the turn-on rate of protogalaxies. Figure 3 suggests 
that i?Q oc {t^f oc (1 + z)-^-^ at z > 2. This gives a decrease by a factor of 2.7 in i?Q from 
z = 3 to 4, consistent with the drop in quasar comoving space density at these redshifts 
(Schmidt, Schneider, & Gunn 1995). 

The three data points that lie above the upper bound in Figure 1 could be explained 
if the enrichment rate /9l = (Fe/H)L/AL of L events were increased by a factor ~ 2. It 
is possible that /3l has a spread of a factor of 2 or 3. Alternatively, we can consider the 
enrichment rate as a steep function of time, with i3\jt) starting quite high and then settling 
down to the value proposed by WQ. In this case (Fe/H) = (Fe/H)p + /.^^ iiiXi')d^ and 
the straight hue evolution in Figure 2a would be replaced by a curve. We note that Fe 
enrichment by Type la supernovae would not be significant for most damped Lya systems 
as such enrichment appears to be significant only at [Fe/H] > — 1 in the Galaxy (Timmes, 
Woosley, & Weaver 1995). 

A more difficult matter is the timescale for condensing protogalactic globs during the 
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expansion of the universe and the timing sequence outhned above. The range in t* required 
is large and indicates that the time required to form protogalaxies and condense most of 
the baryonic matter into stars is comparable to the age of the universe at ^ ~ 2 (~ 3.5 
Gyr after the Big Bang) or possibly even longer. This is in conflict with ab initio models 
that report almost complete condensation of dark matter and possibly cloud or protogalaxy 
formation at 2; ~ 10-20 (e.g., Kamionkowski, Spergel, & Sugiyama 1994). The occurrence 
of [Fe/H] PS —2.6 at 2; ~ 2 implies that there are regions where formation of normal stars 
did not begin until 3.5 Gyr after the Big Bang. The occurrence of [Fe/H] fa —2.6 to —2.4 
in the range of 2; ~ 2.0 to 4.2 indicates that a large fraction of the baryonic matter has not 
been condensed into protogalaxies and stars over the corresponding extended time range. If 
damped Lya systems are random samples of the original medium, this indicates that the 
reservoir of original uncondensed baryonic material may not have been seriously diminished 
over ~ 3.5 Gyr. Indeed, it is possible that the bulk of this baryonic matter is dispersed and 
has not condensed today. Measurements of damped Lya systems at lower redshifts (but at 
times before Type la supernovae contribute Fe) would provide a test. We have no means of 
establishing a priori whether the intrinsic rate of protogalaxy formation decreases with time 
(possibly due to decrease in global density) without appreciably depleting the reservoir of 
baryonic matter. 

A variant of the above scenario is possible which would not be in conflict with the 
models that suggest almost complete condensation of baryonic matter by z ~ 10. The 
modified scenario would be that most baryonic matter in the potential well created by 
non-baryonic dark matter is rapidly collected into Ur-protogalaxies, that formation of very 
massive stars from Big Bang debris is rapid [^2 "C t{z)], and that the explosion of these 
very massive stars destroys the Ur-protogalaxies and redistributes matter into the general 
medium providing a uniform source of Fe (and other elements such as C, N, O, Mg, and 
Si). The material in this second generation medium is then the source for much slower 
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aggregation and formation of protogalaxies and stars. The explosion of very massive stars 
is likely to be very energetic. The potential wells formed by dark matter are considered to 
have a typical escape velocity of ~ 400 km s~^. The velocity of debris ejected from explosion 
of very massive stars is almost certainly much larger than this value. The inter-protogalactic 
medium would then be supplied with hot matter containing Fe, Mg, O, Si, and C. A hint 
in favor of this is found in the recent detection of O VI quasar absorption systems at 
low redshifts (Tripp, Savage, & Jenkins 2000). Knowledge of the explosion dynamics and 
nucleosynthetic products of very massive stars formed from Big Bang debris is fundamental 
to further progress. 

We want to thank Roger Blandford and Marc Kamionkowski for their ongoing interest 
in this work and for generously giving their time and providing, as usual, deep cosmic 
insights and education. Sir Maarten Schmidt has been kind enough to expose us to 
the details of quasar distribution. This work was supported in part by the Department 
of Energy under grant DE-FG02-87ER40328 to Y.-Z. Q. and by NASA under grant 
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Fig. 1. — A summary of the available data on [Fe/H] of damped Lya systems at z ~ 1.5 
to 4.5 (asterisks: Lu et al. 1996, 1997; squares: Prochaska & Wolfe 1999, 2000). The curve 
labeled t"^ = is the upper bound from equation (H) with [Fe/H]p = —2.8, [Fe/H]L = —2.5, 
and Al = 10*^ yr. The curve labeled {t*{z)) = 0.8t{z) appears to represent the average trend 
of the data. Curves for [Fe/H] corresponding to fixed values of t* = 1 and 1.5 Gyr are also 
shown. 

Fig. 2. — (a) Schematic diagram of (Fe/H) as a function of t* at a fixed z. The baseline 
(Fe/H)p is the prompt enrichment value. Normal star formation began in a protogalactic 
system at time after the Big Bang. After reaching (Fe/H)p, a particular mass of matter 
must wait until t* to start the increase in (Fe/H) by contributions from Type II supernova L 
events. For t* = 0, there is a maximum increase in (Fe/H) (trajectory A). As t* approaches 
t{z), the growth is proportionally smaller (trajectories B and C). (b) Schematic histogram 
for the distribution of t* in a set of damped Lya systems at a given z. The distribution for 
z' > z can be obtained by discarding the part of the distribution for z a.t t* > t{z'). 

Fig. 3. — Histogram of t* for data between 2.0 < z < 2.4 in Figure 1 showing a peak at 
t* ^ 2.5 Gyr when t{z = 2.2) ^ 3.2 Gyr. This shows that the rate of turn-on of protogalaxies 
starts from a low value close to the Big Bang and increases until at least z ^ 2. 

Fig. 4.— Values of t* calculated from equation (|) with [Fe/H]p = -2.8, [Fc/HJl = -2.5, 
and Al = 10^ yr for all data in Figure 1. Note the clustering of t* close to t{z) over the 
range of z. Three data points indicated by downward arrows require t* < 0. 
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